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Abstract. - Recently, there has been some vigorous interest in the out-of-equilibrium quasistation- 
ary states (QSSs), with lifetimes diverging with the number N of degrees of freedom, emerging from 
numerical simulations of the ferromagnetic XY Hamiltonian Mean Field (HMF) starting from some 
special initial conditions. Phase transitions have been reported between low energy magnetized 
QSSs and large energy unexpected, antiferromagnetic-like, QSSs with low magnetization. This is- 
sue is addressed here in the Vlasov N — > oo limit. It is argued that the time asymptotic states 
emerging in the Vlasov limit can be related to simple generic time asymptotic forms for the force 
field. The proposed picture unveils the nature of the out-of-equilibrium phase transitions reported 
for the ferromagnetic HMF: This is a bifurcation point connecting an effective integrable Vlasov 
one-particle time-asymptotic dynamics to a partly ergodic one which means a brutal open-up of 
the Vlasov one-particle phase space. Illustration is given by investigating the time-asymptotic 
value of the magnetization at the phase transition, under the assumption of a sufficiently rapid 
time-asymptotic decay of the transient force field. 



Introduction. — Systems interacting via long-range interactions continue to receive a 
considerable interest even in one dimension due to the intricate relationships between their 
dynamical and statistical properties. Of particular interest is the question of the convergence 
of the time-asymptotic dynamics to equilibrium statistical mechanics ensemble predictions. 
This issue is more than purely academic, it is relevant to physical systems ranging from 
hot plasma physics and its ubiquitous wave-particle interactions phenomena, including free 
electron lasers and laser-plasma devices, to self-gravitating stellar systems. The Hamiltonian 
Mean Field (HMF) model [I] 

N 2 N 

ff = Ey + ^E [1-008(9,-9^] (1) 

i— 1 — \ 

is a choice toy model to explore those problems, since it exhibits a non trivial long (and 
actually infinite) range collective dynamics meanwhile permitting a simple exact derivation 
of the equilibrium statistical mechanics. It describes the all-to-all interaction of N parti- 
cles moving on the unit circle with momenta pi and canonically conjugated positions Oi. 
Here c = ±1 contains the information on the attractive or repulsive nature of the interac- 
tion: a positive c gives the attractive ferromagnetic model and a negative c the repulsive 
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antifcrromagnetic one. 

Recently, some unexpected out-of-equilibrium behaviors have been pointed out: For 
the antiferromagnetic non- magnetized case, that is analogous to plasma systems [2J, and 
for low energies, the system spontaneously develops into a biclustered state whose lifetime 
is an increasing function of N [3111]. For the ferromagnetic case, in which equilibrium 
statistical mechanics predicts a second order phase transition, and for some special class 
of initial conditions, the system has been shown to evolve to out-of-thermal-equilibrium 
quasistationnary states (QSSs) with low magnetization, still having lifetimes diverging with 
N, at subcritical energies that would be thermodynamically associated to the magnetized 
phase. This means that, depending on initial conditions, observables in the Vlasov N — > oo 
limit of the V-dimensional dynamics ([TJ do not converge asymptotically in time towards 
their equilibrium statistical mechanics predictions: the limits N — > oo and t — > oo do not 
commute [5HI]. 

The Vlasov equation that governs the evolution of the distribution function f(p,9,t) 
forms then the natural framework to investigate those QSSs. For the HMF model in the 
ferromagnetic case (c = 1), it reads 

-+p- + E(9,t)-=0, (2) 

where E stands for the force field given by 

+ 00 2-7T 

E(9,t) = - f du J dasm(9-a)f{u,a,t). (3) 

-oo 

Eq. ([3]) forms a consistency relation that makes Vlasov equation @ intrinsically nonlinear. 
It is meaningful to introduce the magnetization vector M = (M x ,M y ) through M x [/] = 
/ / dudacosaf (it, a, t) and M y [f] = j f dudasmaf (u,a,t), which enables to write the 
one-particle Hamiltonian associated to ([2j) as 

fn (p, 0, i) = |- - M x [/] (i) cos 9 - M y [/] (t) sin 9. (4) 

The first aim of this Letter is to propose that the time asymptotic states emerging 
in the Vlasov limit be related to simple generic time asymptotic forms for the force field 
E (9,t). The proposed picture unveils the nature of the out-of-equilibrium phase transitions 
reported for the ferromagnetic HMF: This is a bifurcation point connecting an integrable 
Vlasov time-asymptotic one-particle dynamics to a partly non-integrable and ergodic one. 
Illustration is given by investigating the time-asymptotic value of the magnetization at the 
phase transition, under the assumption of a sufficiently rapid time- asymptotic decay of the 
transient force field. Then, a proof of principle is given, that the time-asymptotic state 
of the system can be completely derived as an initial value problem, in the domain where 
initial magnetizations are sufficiently large that the out-of-equilibrium phase transitions be 
of second order. This serves to specify the validity regime of the proposed time-asymptotic 
picture. 

Numerical evidence. — Considering the ferromagnetic Hamiltonian (fTJ) with c = 1, it 
can be shown [TJ that the system undergoes a second order phase transition for U = H/N — 
e c = 3/4 with order parameter the modulus of the magnetization M = ||M||. QSSs have 
been observed starting with initial waterbag distributions of particle momenta and positions 
of the form 

h(p,8) = 4 . ^ 1 [- Ap; Ap] 0) 1 [- A0; A9] (<?) , (5) 

where tx denotes the characteristic function of the domain X, Ap > and < A9 < it. 
There is then a one-to-one relationship between (Ap, A9) and (U,Mq), where Mo is the 
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initial magnetization, through M = sin (A0) /A0 and U = Ap 2 /6 + (1 - M§)/2. To be 
specific, a transition between two sorts of time-asymptotic states has been reported. For 
a given A0, that is for a given initial magnetization M , there exists some critical energy 
density U c , whose value is a growing function of M , such that: For U < U C (M ), the 
system converges asymptotically towards the clustered state as predicted by equilibrium 
statistical mechanics whereas, for U > U c (Mq), magnetization drops and the system mimics 
an antiferromagnetic (upper critical) behavior in the limit of large N . This latter regime 
displays a specific signature: whereas, for U < U c (Mq), the one-particle phase space is that 
of the usual single resonance pendulum centered on p = associated to the clustered phase, 
for U > U c (Mq), it shows a superposition of two waves traveling at opposite phase velocities. 
This has been enlightened very recently through very detailed numerical simulations in Ref. 
[8] (see also Ref. [9] for large time one-particle phase space plots obtained from Vlasov 
simulations of the HMF model in the low- magnetized phase). 

These out-of-equilibrium phase transitions have been classified in two types: for a suffi- 
ciently low initial magnetization, below some threshold value Mot ~ 0.15, transitions are of 
the first order type with the order parameter, namely the magnetization, being discontin- 
uous at the critical point. For a larger initial magnetization Mq > Mot, transitions are of 
second order with no discontinuity of the magnetization at the critical point. This has been 
reported in extensive finite- N numerical simulations in Ref. |10j and the agreement with 
statistical predictions by P.H. Chavanis [IT] following Lynden-Bell's approach has proved 
to be remarkably good. However, a careful observation of QSS magnetizations in Fig. 2 
shows that, at least in the second order regime of out-of-equilibrium phase transitions for 
Mo > Mot (bottom plot), its upper critical values appear to remain strictly positive even 
in the N — > oo limit. This means that there is a continuity of the order parameter at the 
transition point yet this does not vanish in the high energy phase. This is absolutely consis- 
tent with the above mentioned QSS phase space pictures. Actually, the observed long time, 
upper critical, phase space inhomogeneity reflected by the two contra-propagating clusters 
implies that the modulus of the magnetization is, even small, strictly positive and that 
particles do not behave asymptotically in time as free particles even in the Vlasov limit. 

This latter behavior belongs to a more general phenomenology. It is similar to the one 
emerging commonly from long time simulations of antiferromagnetic-like ID Vlasov-Poisson 
systems [13]. Unless the electric field is totally damped to zero, there is indeed a shared 
agreement that time-asymptotic states can be represented as a superposition of traveling 
waves [T5[I16|. The physical picture behind this is that nonlinearities, that eventually come 
into play, take place in the form of particle trapping which freezes the dynamics. 

Time-asymptotic bifurcation. — By analogy with the almost periodic time-asymptotic| 

states observed in plasma Vlasov-Poisson systems, it is possible to infer that the time asymp- 
totic form of the force field, A, may be put in the form of a sum of traveling waves. Having 
this in mind, the time asymptotic force field corresponding to the usual magnetized clustered 
phase, with U < U c (Mq), may be written as a time-independent (zero phase velocity) single 
wave 

A{6) = -a sin (0 - tp) . (6) 

On the contrary, for U > U c (Mq), the system displays a typical plasma-like (i.e. antiferro- 
magnetic) behavior and the time asymptotic force field may be written as a superposition 
of two traveling waves with opposite phase velocities u> > and — u 

A (0, t) = sin (0 - p + <p - Lot) - - sin (0 - p - <p> + uit) . (7) 

It can be simply checked |12) that this two-wave form (JJJ is the minimal superposition of 
traveling waves compatible with the HMF cosine potential ([3]). Here ip and <p represent 
constant angles and a > is the constant wave amplitude. All these parameters (a, w, p 
and p) relate to the initial conditions through U and M , yet this dependence has not been 
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Fig. 1: Schematic view of the out-of-equilibrium bifurcation in terms of wave pulsation between 
asymptotic ferro- and antiferromagnetic states. Only the uj > half-space has been represented. 



made explicit to ease notations. In order to continuously match (j6|) and at the transition 
point, we should take <p = and u> — + in when U — U c (Mq). Eqs. ([6]) and (0 capture, 
so to speak, the asymptotic skeleton of the Vlasov HMF dynamics evolving from the initial 
conditions (0. 

The continuity of u> at the transition is an assumption compatible with available nu- 
merical simulations undertaken in the second order regime of the out-of-equilibrium phase 
transitions. For instance, the bifurcation diagram plotted in Fig. 3 of Ref. [8] shows that 
the phase transition is signaled by a pinch in the velocity space. As the energy increases as 
one goes right from this pinch, two resonances first fully overlap then gradually emerge in a 
symmetric way. This is further commented in the concluding discussion. A scheme of this 
time-asymptotic bifurcation is depicted in Fig. [TJ The out-of-equilibrium phase transition 
is then associated to the adiabatic limit w — > in Eq. ([7]) that separates an integrable one- 
particle dynamics for the condensed phase to a one-and-a-half degrees of freedom (d.o.f.) 
one-particle dynamics for the antiferromagnetic-like phase |14j . 

Analysis. — Investigating the time asymptotic fate of the system as an initial value 
problem forms a formidable task, because of the nonlinearity of the Vlasov equation induced 
by ((3|), that is out of the scope of the present paper. Instead, we wish to investigate further 
the implications of the picture presented in Fig. [T] to unveil the dynamical nature of the 
out-of-equilibrium phase transition. To do so, let us follow briefly the approach proposed 
by Lancellotti and Doming (TSHU] for the ID Vlasov-Maxwell system, and decompose the 
force field E into its time-asymptotic part A and a transient part T such that 



with lim t _ i . 00 T {9, t) — 0. Our interest lies, in particular, in the time asymptotic value of the 
magnetization, that is a collective observable, in the vicinity of the phase transition. Due 
to the decomposition (j8]) it is possible to extract a in Eqs. ([6]) and ([TJ from suitable time 
averages. This can be straightforwardly done through 



E(9,t) = A(6,t)+T(6,t) , 



(8) 



t 2-rr 
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that gives, using ([3]) and ([8]);, 

l + F 7 \ 

2 (1 + <^w,o) = ( du daf (u, a, r) cos (<^5 — wt) cos (a — ip) ) , (10) 

Loo ' r 

where (.) T = lim^oo \ J Q dr denotes the time average. In Eqs. © and (fTUI) . we used 
(p{u) = 0) = 0. As previously mentioned, this, together with the continuity of ui at the 
transition, can be called a "second order type" assumption in the sense that one moves 
continuously from the low energy to the high energy phases at the critical point. 

In order to relate the time- asymptotic characterization to the initial state, a rather crude 
assumption will be done since the contribution coming from the transient force-field will be 
discarded. Practically, this amounts to consider that the decay of the transient force-field 
is so fast that the form of the particle distribution function does not deviate substantially 
from that of the initial condition, so that one can write for any (p, 9, i) 

I ( P , e, t) ~ f [ P i ( P , e, t) , e£ ( P , e, t)] , (11) 

where fo {p 1 0) = f (p, 9, t = 0) stands for the initial distribution function and where [pf (p, 9, t) , 9f (p, 9, £)]| 
denotes the phase space location at time s of a particle arriving in (p, 9) at time t un- 
der the time-asymptotic one-particle dynamics. This derives from the one-particle time- 
asymptotic Hamiltonian h\(p,9,t) through the characteristics d9/dt = p = dh\/dp and 
dp/dt = A(9,t) = -dh!/d9. 

Nature of the out-of-equilibrium phase transition. From now on, we focus on 
the lo — limit that is associated to the out-of-equilibrium phase transition and look for an 
illustration of the implication of its dynamical signature as depicted on Fig. [T] Due to Eq. 
(fTTj) . the identity (flU)) becomes simply 

du J da fo [uq (u, a, r) , (u, a, r)] cos (a — ip) J . (12) 
v -oo o ' T 

Using Liouville's theorem on the conservation of phase space and Eq. ([5]) and inverting the 
order of the time integration, this reads 

+Ap A0 

1 



4ApA6» 



du j da (cos (a - ip)) T , (13) 

Ap -A8 



where the index A has been dropped to ease notations. 

Limit uj — > 0~. Let us first evaluate (fT3|) from the left (integrable) side. Then, since 
A does not depend on time, the one-particle Hamiltonian hi(p, 9) = p 2 /2 — a cos (9 — ip) is 
integrable, and the time average in (I12|) can be replaced here by an ensemble average on 
the energy level h\(u,a) = hi(uo,ao) = Eq. Let us denote by Q(x) the usual Heaviside 
function defined by 0(x) = J_ dy5(y) and by Z(a,E) the one-particle volume partition 
function defined by 

OO 7T 

Z(a,E)= J du J daQ[E-hi(u,a)]. (14) 

— OO — 7T 

Then, it is easy to check that 

1 f> 7. 

(15) 



(cos(a-<p)) T = I d " Z 
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For the trapped motion —a < Eq < a, moving to action-angle variables yields Z(a, E) — 
J J dJdip®\E-h x {J)\ = 2tt / dhuj(h)~ 1 Q (E — h) with u(h) = Ky/a/(2K(k)) [E] where 

k = y / (h + a) /2a and K(k) = J*' 2 d<p/ \J\ — k sin 2 0, the complete elliptic integral of the 
first kind. This gives, using (TT5| . 



2E( 



Eo+a-] 
2a 



(cos(a-y>)) r = „} E ™ a ( -1, (16) 

2a 



where E{k) = f* /2 d^l - k sin 2 <p is the complete elliptic integral of the second kind. For 
the untrapped motion Eq > a, one has directly Z(a,E) — 2y/2 da^J Eq + a cos (a — ip). 
This yields 

/ I w ^ + 1 ) J? (^) E Q 

(cos (a - p)) T = . \ ^ . (17) 



( 2a 
\a+_E, 



^ ' i+E a ) 

It remains to evaluate the double integral (TIB"]) using the expressions (TT6"|) for — a < E < a 
and pT|) for B > a, with _E = u^/2 — a cos (a — ip). 

Limit w —> + . Let us now consider the limit lu —> + of the 1.5 degrees of freedom 
one-particle Hamiltonian with force field given by Eq. ([7]). This can be written as A(9, t) = 
— a cos (ip — cot) sin (9 — ip). This corresponds to the pendulum with variable amplitude. 
The essential thing to recognize here is that the limit to — > is singular: Actually, because 
of separatrix crossings, the 1.5 d.o.f. Hamiltonian associated to ((?]) exhibits an ergodic, 
diffusive dynamics on a confined phase space region even in the adiabatic limit uj — > + . 
The ergodic region corresponds to the phase space region swept by the resonance cat's eye 
in which particles transit between trapped and untrapped motion [18H20) . Consequently, 
in the limit w — > + , the evaluation of a in Eq. (|13p differs from the previous lu — > CP 
calculation because of a different contribution from the initial phase space domain swept by 
the separatrix, i.e. such that —a < hi(uo,ao) < a. In this ergodic domain, (cos (a — tp)) T is 
given by its average on the inner whole cat's eye region 



/ da cos (a — ip) y/l + cos (a — ip) 

(cos(a-p)) T = - ^ = ^— = ~. (18) 

J da.yf \ + cos (a ~ p) 

— TV 

The outer region, for which hi(uQ,oto) > a, corresponds to a regular one-particle motion 
and gives the same contribution to (|13[) than in the lo — > 0~ case. In evaluating the double 
integral (TT5]) in the lu — >• + limit, one should then use the identity pT|) for Eq > a, while 
using the expression (TT8l) for — a < Eq < a. 

Implications. The lu — > 0~ and lu — > + values of the time-asymptotic magnetization 
a were numerically computed. Figure [2] represents the curves, in the (Ap/s/a,ip) space, on 
which the time asymptotic magnetizations a for lu = 0~ and lu = + coincide, for several 
values of the initial magnetization Mq. Here the range of the angle <p is restricted to the 
interval [0;tt]. The full plot in the [0; 2ii\ range can be recovered by an axial symmetry 
about the horizontal line tp = n (since when tp is solution, then 2ir — <p is also solution by 
switching the up and down waves in (J7J). When Ap/^/a becomes larger than 2, <p remains 
constant since resonance cat's eye has been fully covered and that regular contributions (|17[) 
from the outer region cancel out. 

Interestingly enough, when Mq is sufficiently small (roughly below 0.15), there appears 
a gap, namely a region of the rescaled parameter Ap/y/a, where the time asymptotic mag- 
netizations for lu = and lu — + do not intersect. This is reminiscent of the first order 
nature of the out-of-equilibrium phase transitions, associated with a discontinuity of the 
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Fig. 2: Curves in the (Ap/i/a, tp) space on which a(ui — > ) = a(uj — > + ) labeled by the values of 
the initial magnetization Mq. The range of the angle <p is restricted to the interval [0; tt]. 



order parameter at the transition, evidenced in some already mentioned previous numerical 
works 10; and thermodynamical analysis [IT] , for the same domain of Mq. It is interesting 
to note here that the analysis by itself is able to unveil that the presumed continuity between 
Eqs. ((6])-([TJ) cannot be satisfied in that domain, since a is not continuous. 

Derivation of the time-asymptotic state: a proof of principle. In order to 
fully characterize the transition, namely, for a given Mq, determine the time-asymptotic val- 
ues of the magnetization a, of the wave phase ip and determine the transition energy through 
Ap, an extra condition is needed. The natural additional relation to fulfill is given by the 
continuity of the energy at the transition. Its expression is h[f) = K[f] + 1/2 1 — (M[/]) 2 , 

with the kinetic energy K[f] = J J p 2 /2/ (p,8,t) dpdd. In the regular domain, namely for 
the case ui — > 0~ with h\{p, 9) = hi (p , 9 ) = p 2 /2 — a cos (9 — ip), neglecting the effect of 
transient and using Liouville theorem, the time average of the kinetic energy is given by 

+Ap AS +Ap AO 

= ZA^Atf / dpo J d6ohl{pM+ A^Ae J dp ° J ^Mcosfl? -*>)>« 

-Ap -A8 -Ap -Ad 

+Ap A6 

^ 2^ P J dp 4 ~m I ^cos(^) + a 2 , (19) 

-Ap -A6 

where the identity (|13[) was used. The continuity of the energy at the phase transition 
amounts then to the identity 

£f - aMocos^ + a 2 - (K[f])r° + = \ (<M[/] 2 >^°" - (Af[/] 2 )r° + ) , (20) 

where time averages may be computed using the exact characteristics for h\ in the u> — > _ 
case and the numerically obtained ones in the ui — > + case. 

This procedure would amount to determine the transition as an initial value problem, a 
truly desirable perspective that would involve the knowledge of the sole one-particle Vlasov 
dynamics without any reference to equilibrium statistical mechanics. Yet this desirable 
perspective is out of reach of the present analysis, surely due to the presently too stringent 
condition on the decay of the transient force field. The obtained values of the transition 
energies |12j are actually sensibly below the numerically [8] or thermodynamically computed 
ones [IT]. This discrepancy can be already inferred from Fig. [3] This Figure displays the 
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possible time-asymptotic values of the magnetization a at the phase transition in the case 
where a(u) 0~) and a{uj — > + ) coincide. As in Fig. [3J it is clear that, for low enough 
M , there exists some forbidden range for Ap that may signal a discontinuity of the order 
parameter associated to the first order transition region. It is also clear from this picture, 
that, at least for low enough Mq, the possible Ap are bounded by values (~ 1) which would 
yield values of the energy transition U C (M ) below the available numerically computed [10] 
and thermodynamically predicted ones . Finally, in the domain of initial magnetizations 



a 




0.2 0.4 0.6 0.8 



Fig. 3: Curves in the (Ap, a) space on which a(cj — > ) = a(uj — > + ) labeled by the values of the 
initial magnetization Mo. 

where the out-of-equilibrium phase transitions are of first order, the additional identity ([20| 
is insufficient to close the system of time- asymptotic parameters because a is not continuous 
at the transition. This signals that the continuity assumption made in Eqs. © and (|7|) 
breaks: uj(U, Mq) is strictly positive in the upper critical domain and its knowledge at the 
" right side" of transition line is required. 

Conclusion. — To the author's knowledge, this study is the first to address the issue 
of out-of-equilibrium phase transitions in the HMF model from the perspective of an initial 
value problem in the Vlasov limit. The essential point of this Letter is to propose a picture 
unveiling the nature of the phase transition: This corresponds to a brutal open-up of the 
time-asymptotic Vlasov one-particle phase space. The phase transition coincides with a 
jump between a time-asymptotic phase space foliated by energy lines and a phase space 
divided between an ergodic and a regular [21] components. It is proposed as a mechanism 
for second order phase transitions compatible with non-vanishing time-asymptotic values 
of the order parameter in mean-field long-range systems. As a byproduct, this study shows 
indirectly that the a priori reasonable hypothesis of a rapid decay of the time-asymptotic 
transient force field may not be satisfied by the HMF model. This is a posteriori not 
astonishing in view of the numerous recent papers reporting incomplete relaxation and 
deficient mixing properties in such a long-range system [22j - !26| . It remains to elucidate 
whether relaxing the hypothesis on the transient field would still be compatible with a purely 
dynamical initial value treatment, e.g. in the spirit of Lancellotti and Dorning's approach. 
The present approach underlines the importance of the knowledge of the time-asymptotic 
spectrum w{U, Mq) of the dynamics. It would be instructive to proceed to time-Fourier 
transforms of large time Vlasov numerical simulations in order to extract this information. 

Very recently, the author became aware of the just published paper by Leoncini et al. 
[27] whose possible connections with the present approach would be interesting to explore. 
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